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Generalized qubits of the vibrational motion of a trapped ion 
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We present a method to generate qubits of the vibrational motion of an ion. The method is 
developed in the non-rotating wave approximation regime, therefore we consider regimes where 
the dynamics has not been studied. Because the solutions are valid for a more extended range of 
parameters we call them generalized qubits. 
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Nonclassical states of the center-of-mass motion of a trapped ion have played an important role because of the 
potential practical applications such as precision spectroscopy [1] quantum computation [2, 3] and because of funda- 
mental problems in quantum mechanics. Ways of producing Scrhodinger cat states [4], squeezed states [5], nonlinear 
coherent states [6, 7], number states, specific superpositions of them, and in particular, robust to noise (spontaneous 
cmmision) qubits have been proposed [8] . In theoretical and experimental studies of a laser interacting with a single 
trapped ion it has been usually considered the case in which it may be modeled as a Jaynes-Cummings interaction 
[5, 9, 10], then exhibiting the peculiar features of this model like collapses and revivals [11], and the generation of 
nonclassical states common to such a model or (multi-photon) generalizations of it [12-14]. In treating this system 
I . usually two rotating wave approximations (RWA's) are done (the first related to the laser [optical] frequency and 
the second to the vibrational frequency of the ion), to remove counter- propagating terms of the Hamiltonian (that 
can not be treated analytically). Approximations on the Lamb-Dickke parameter, 77, are usually done, considering 
it much smaller than unity. Additionaly, other approximations are done, based on the intensity of the laser shining 
on the trapped ion: the low-excitation regime $7 <C v and the strong-excitation regime 51 3> v [15], with il being the 
intensity of the field, and v the vibrational frequency of the ion. 
-H Recently there has been an alternative approach to the study of this dynamics [16]. In this approach a unitary 

J> ■ transformation is used in order to linearise the ion-laser Hamiltonian. This transformation has been also used to pro- 
pose schemes for realising faster logic gates for quantum information processing [3] . Under this unitary transformation 
the Hamiltonian takes exactly the form (note that not an effective form) of the Jaynes-Cummings Hamiltonian plus 
an extra term (an atomic driving term). In such a case a RWA may be done [16] in order to have an analytical solution 
for this problem, but it brings with it a new condition: 17 of the order of 1v (note however that this is a new regime) 
and r\ still much less than unity. Later, another transformation [17] was used to diagonalize the linearized ion-laser 
Hamiltonian, without further conditions on 17 or r\. This allowed the diagonalization of the Hamiltonian only in the 
ion basis. Exact eigenstates of the ion-laser Hamiltonian, i.e. trapping states for this system have been found [18], 
but because they do not form a basis, a complete (exact) solution may be found only for such states (eigenstates). 
In this contribution we consider the complete Hamiltonian for the ion-laser interaction, linearise it as in [16] and 
k^ ' further unitarily transform it, without performing the RWA to obtain an effective Hamiltonian that can be easily 
solved . This is an extension of a method of small rotations recently applied by Klimov and Sanchez-Soto [19] to the 
Cy_' Dicke model and other systems (they apply small rotations on the atomic basis). This allows us to obtain a more 
general solution valid in a more extended range of parameters. We apply this solution to a simple initial state and 
show that qubits may be produced. Because may be produced with less constrains on the parameters, we call them 
generalized qubits. 

We consider the Hamiltonian that describes the interaction of a single two-level trapped ion with a laser beam (we 
set h = 1) [9, 15] 

H = im + -»„ + £7(<T_e-"'( a+at > + a+e 1 ^^), (1) 

where h — a^a, with a* ( a ) the ion vibratonial creation (anihilation) operator, and <7 + = \e)(g\ (<r_ = |<7)(e| ) is 
the electronic raising (lowering) operator, |e) (\g)) denoting the excited (ground) state of the ion. The detuning S is 
defined as the difference between the atomic transition frequency (luq) and the frequency of laser (wl). 
Applying the unitary transformation T\ [16] 

f x = -L Q[£t(/9) + b{fi)]I + i[Dt(/9) - D((3)}a z - &(/3)<r_ + D((3)a + ^j , (2) 



(where I = \g)(g\ + |e)(e| and D(f3) is the displacement operator with j3 = ^ its amplitude) to (1), we obtain the 
linearised Hamiltonian (in the on resonance case i. e. 5 = 0), given by [16] 
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Hx = TiHfl =uh+ Vta z + i^-(a - a 1 )^- + &+) + ^-. (3) 

We now apply the unitary tranformation 

j" 2 = e -ie(a+a t )(CT + +ff_) /^\ 

to (3) to obtain H 2 = f 2 H\fl, 

H 2 = v(h + ie{a-a)){d-+d + ) + e 2 ) 

+Vt (a z cos[2e(a + a f )] + z(er_ - <7+) sin[2e(a + a))]) 

+i^ ( a-at) ( a_ + a + ) + ^. (5) 

By considering e <C 1, we can approximate (5) as (we disregard constant terms that only contribute to a shift of 
energies) 

H 2 ~ v(h + it{a-a r ){a- +<7+)) + £l(a z + 2ie(a- -o' + )(a + a t )) 

+i^(a-at)(a_+<7 + ), (6) 



and by setting 



we finally obtain 



C 2*/ + 2fi' lJ 



i?2 = w + £7<r 2 + zA(fi + a — a^CT_). (8) 



Note that the coupling constant, A = JT^n > nas changed (before it was ^). 

We should remark that transformation (4) requires e< 1 and this may be achieved in different forms: 

a) r\ <C 1 and v and f2 any numbers (this is, it can be v <C f2, f ^> 17 or of the same order of magnitude); 

b) no restrictions on r\ and ^ ^C 17 or, 

c) 7] < 1 and f < 17 (note that, for instance, a value of n = 0.3 and 17 = 2v gives e = —0.03). 

The three possibilities above allow the Hamiltonian to be approximated to first order and to disregard terms of 
second and higher orders. 

We are now in the position to give a solution to the Hamiltonian (1), that reads 

|*(*)> = f ^UT\^(0)), (9) 

where we have written 



T = T 2 T Y 

1 (I 

~ 71 



Q[£t08_) + b{P-)]I + \\&<jS-) D((3_)]a z Dt(/3_)a_ + D(/?_)cr + ) , (10) 



with j3- = i(% — e) and where 1^(0)) is the initial wave function and U is the evolution operator of the off-resonant 
Jaynes-Cummings Model 

fj _ „—it(vn+\va z ) -it[Aa z +i\(aa + -a f a-)] Qj\ 

where we defined A = £7 — %. Equation (11) may be finally written in the form [20] 

U = e-^^+i"*') (^[Un + U 22 ]I+ i[J7n - f7 22 ]<r z + *7 21 a_ + ?7 12( t + ") , (12) 



with 



and 



fr - ± . A sind fi+ it 

Uii = cosaft + ir — ill — , (16) 

U 12 = \a — , (14) 

aft 

7T t sind A+ it 

C/ 2 i = -Aa 1 — : , (15) 

Oin+l 



smctnt 

U22 = cos afit + z A — , (16) 

a?,. 



where 

d ft = vam"^. (17) 

By considering the ion initially in its excited state and in a coherent (vibrational) state with amplitude /3_, i.e. 

|tt(0)> = |e)|/3_), (18) 

we obtain the evolved wave function at time t 

\Mt)) = \t){BJ) ( [e'^ (cos cat -i— sinc^t) + e iO *]|0) + — e _< ^ sinaitll' 
2 y aj ai 



+ -£» f (/3-) I [e- 4 ^(cosait-i— sinaii)-e tOt ]|0)-— e-^sina^ll) ) | ff ) 
2 \ ai ot\ J 



(19) 



with ai = (l|dft|l). Note that D(/3_)|fc) = \/3-,k) is a displaced number state [21]. By measuring either the 
excited state or ground state of the ion (there are standard techniques to do so, see for instance [9]), we end up 
with a superposition of displaced number states with amplitude /?_ or — /3_, respectively. Therefore, by displacing 
the resulting state if the ion is measured in the excited state by — /3_ (or by /?_ if measured in the ground state, we 
produce a qubit of the vibrational motion of an ion. To clarify this point, let us assume the ion is measured in the 
excited state. The vibrational wave function collapses to the state 

\®vib(t)) = -^D(P-) ([e-'Tr(cosai£-i— sinarf) + e mt }\0) + — e _i ^ sinaitll) ) , (20) 

N \ a\ a\ J 

where N is the normalization constant. By finally displacing the state by an amplitude — /?_, we obtain the (qubit) 
state 

|* d (t)> = — I \e- li £ (cos arf-i— sinaii)+e jm ]|0) + — e _i ^ sinaiill) ) . (21) 

N \ ot\ a.\ J 

Moreover, without (conditional) measurement, and by controling the interaction time between the ion and the laser 
beam, setting it to ait = 7r, the (Schrodinger cat) state 

Mait = tt)) = I ((e int - e-^)\^)\e) - (e int + e~ ivt )\ - 0_)\g)) , (22) 

is generated. This state was produced in [4] and also studied in [22]. Note that the quantity j3- is slightly larger than 
the Lamb-Dicke parameter, which in general is not needed to be small. 

In this contribution we have suceeded in solving the problem of trapped ion interacting resonantly (6 = case) 
with a laser field whitout making the rotating wave aproximation, but applying a method of small rotations that 
allowed such solution. It should be remarked that now there are new regimes that can be reached, since the usual 
approximations as the low intensity or high intensity regimes were not considered. 
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